Friction, order, and transverse pinning of a two-dimensional elastic lattice 
under periodic and impurity potentials 



O 

o 
o 

(N 



J2 

o 

CD 



Takaaki Kawaguchi 
Department of Technology, Faculty of Education, Shimane University 
1060 Nishikawatsu, Matsue 690-8504, Japan 

Hiroshi Matsukawa 
Department of Physics, Osaka University 
1-16 Machikaneyama Toyonaka, Osaka 560-0043, Japan 

Frictional phenomena of two-dimensional elastic lattices are studied numerically based on a two- 
dimensional Frenkel-Kontorova model with impurities. It is shown that impurities can assist the 
depinning. We also investigate anisotropic ordering and transverse pinning effects of sliding lattices, 
which are characteristic of the moving Bragg glass state and/or transverse glass state. Peculiar 
velocity dependence of the transverse pinning is observed in the presence of both periodic and 
random potentials and discussed in the relation with growing order and discommensurate structures. 
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Static and dynamic properties of an elastic or plastic 
media, subject to an external force and under periodic 
or random potentials, have attracted much attention re- 
cently in relation to the vortex lattice in type-II supercon- 
ductors, charge and spin density waves, Wigner crystals, 
the interface growth, the dry friction between two solids 
and so on They are pinned by those potentials 

under a small external force and start a kind of sliding 
motion above a certain threshold force. Various peculiar 
phenomena are predicted and/or observed in these sys- 
tems. Most theoretical or numerical works study, how- 
ever, systems under only random potentials or those un- 
der only periodic potentials. The latter ones are investi- 
gated as models of vortex lattice with intrinsic pinning, 
commensurate density waves, the friction between atomi- 
cally flat surfaces, etc. In actual experiments of those sys- 
tems, a finite amount of randomness will exist. Theoret- 
ical or numerical works of systems under random poten- 
tials in addition to periodic potentials are indispensable. 
The behavior of friction in the presence of impurities in 
addition to the periodic potential has been studied, based 
on a simplified one-dimensional model, numerically [^]]|] 
and theoretically |ic|| . The interfaces where the friction 
acts are, however, in most cases two-dimensional, and 
the frictional properties are expected to strongly depend 
on the dimensionality of the system. In this paper, we 
investigate numerically the behavior of friction based on 
a two-dimensional elastic lattice in the presence of im- 
purity and incommensurate periodic potentials ]ll| ] . The 
model is also applicable to the vortex lattice in thin films, 
and so on. We show that the maximum static frictional 
force in the presence of impurities can be lower than that 
in the absence of them. This is caused by the frustration 
between impurities and the periodic potential. We also 
investigate the transverse pinning [^|-|5). When we apply 
an external force along the transverse direction to the 
sliding lattice, the sliding velocity has a finite transverse 



component only for the transverse external force that is 
greater than a finite critical value. This transverse pin- 
ning is predicted in the case under random potentials 
alone and the manifestation of a quasi long-range order 
(QLRO) along the transverse direction. We calculated 
the structure factor of the sliding lattice, which show 
clearly such transverse QLRO. The velocity dependence 
of the transverse threshold force can show nonmonotonic 
behavior in the presence of both periodic and random 
potentials. The velocity dependence of the kinetic fric- 
tional force is also presented. They are discussed in the 
relation with the QLRO and discommensurate structures 
of sliding lattice. 

As a simple model which describes an interface be- 
tween two solids, we consider here a two-dimensional 
Frenkel-Kontorova (FK) model with impurities. 

The model consists of a square lattice of atoms in- 
teracting with each other via a harmonic force on 
a substrate. The periodic potential resulting from 
the substrate is represented by a two-dimensional si- 
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nusoidal form, Up(r) 

The impurity potential t/j(r) is expressed as, Ui(r) 

~ "T J2i 1 ex P (~4(r — Rj) 2 ); where Nj is the number of 
impurities and the coordinate of the ith impurity dis- 
tributed randomly on the substrate. The overdampcd 
equation of motion of each atom is given by 



+F P <(u iij )+F?(u i , j ) + 
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where Ujj is the position of the atom with index numbers 
i and j, and the superscript a stands for x and y. In 
the above, F P {r) = -VU P (r) and Fj(r) = -Wz(r) 
are the force from the periodic potential and that from 
impurities, respectively, and F ex is the external force. 
We apply F ex along the x-axis, except in the study of 
the transverse pinning. 
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In this paper we consider an incommensurate system 
and choose the ratio between the mean atomic spacing 
of the lattice, c a , and the period of the periodic poten- 
tial, c p , as the golden mean ratio: c a /c p = \/5 + 1/2 = 
1.618 • • • (we set c p = 3 and c a — 1.618 • • • x c p ). The 
number of atoms in the samples N is 55 2 , 89 2 , or 144 2 , 
and the number of impurities iVj is set so that the impu- 
rity density Ni/(Nc 2 ) is fixed at about 0.21. Numerical 
simulations are performed using a Runge-Kutta formula 
under periodic boundary conditions. Hereafter, v a is the 
velocity averaged over space and time. 

Note that the periodic boundary conditions force the 
system as a whole to remain incommensurate. The re- 
sults obtained here are, however, expected to be semi- 
quantitatively the same for a large system with open 
boundary conditions which is originally incommensurate, 
such as experimental systems. A kind of discommen- 
surate structure with commensurate domains separated 
by defects plays an essential role for the present results. 
Such a structure appears also for the large system with 
open boundary conditions. Hence we can simulate the 
large system with open boundary conditions by a small 
system with periodic boundary conditions. 

First, we investigate pinning properties of static states 
(v x = v y = 0) by evaluating the depinning threshold 
force F c for various values of K and W. According to 
the scaling theory of the pinning of density- waves p3[ , 
F c oc W^ 1 for weak pinning and F c oc W for strong 
pinning, where d is the spatial dimension [ljj ]. On the 
other hand, in the presence of only the periodic poten- 
tial, the lattice is pinned due to the occurrence of the 
Aubry transition |H for K > K c , where K c w 0.23 in 
the present system. We can expect that these two pin- 
ning mechanisms compete with each other, as observed 
in one-dimension |p|-]l0|. 

In Fig. 1 the depinning threshold forces F c 's are plot- 
ted against W. In the absence of the periodic potential, 
K = 0, the weak pinning behavior, F c cx W 2 , is clearly 
observed in a small W regime, W <C 10, and then, with 
increasing W, the F c — W curve gradually approaches the 
strong pinning behavior, F c oc W. When the periodic 
potential with a finite K (> K c ) is present in addition to 
impurities, however, the weak pinning behavior is mod- 
ified especially in the small W regime. In this case, the 
pinning is governed dominantly by the periodic potential 
and the effect of impurities becomes only a weak pertur- 
bation to the pinned state. As the strength of impurity 
potential W is increased, the pinned state crosses over 
to the weak or strong pinning regime. In the crossover 
regime, an interesting behavior is observed, that is, F c 
takes a minimum value. In this regime, impurities dimin- 
ish the pinning effect of the periodic potential. In other 
words, impurities assist the depinning. Such a behavior 
is observed in a one-dimensional system || and caused 
by the frustration between impurities and the periodic 



potential fX0|. 

Next we focus on sliding states. We consider here 
cases in which the impurity potential is much stronger 
than the periodic potential, W ^> K, for which the 
static states are pinned dominantly by the impurity 
potential. To clarify the lattice order, we calculate 
the structure factor in steady sliding states: S(k) = 
j expik • (iij — Uj). Typical results of S(k) at 
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some sliding velocities are shown in Fig. 2. In a low 
velocity regime (Fig. 2(a)) Bragg peaks at k x ^ dis- 
appear. This state corresponds to the moving transverse 
glass (MTG) state discussed in Refs. where QLRO 

only along the transverse direction exists. So far, the 
MTG states have been observed in numerical simulations 
of vortex lattices, where the longitudinal positional long- 
range order is destroyed by the dislocations with Burger's 
vectors parallel to the moving direction Q. Such dislo- 
cations are produced dynamically by the random force 
resulting from the interaction with impurities. Also, for 
the present harmonic lattice case such disordered sliding 
states appear. Even in small W cases, when the velocity 
is sufficiently small, the MTG or a MTG-like anisotropic 
state can be observed. The appearance of the MTG state 
would be a common feature for two-dimensional driven 
lattices under random potentials. As the driving force 
is increased, Bragg peaks at k x ^ appear and grow 
gradually, i.e., the MTG state crosses over to the moving 
Bragg glass(MBG) state, where QLRO exists along both 
directions. In Fig. 2(b), the satellite Bragg peaks on 
the k x axis (k x ^ and k y — 0) are slightly higher than 
those on the k y axis (k x = and k y ^ 0). The width 
of the peaks is also anisotropic and broadens in the k v 
direction. The anisotropy of the MBG survives in the 
high velocity regime. Note here that, in two-dimensional 
systems, the change from MTG to MBG would occur as 
a crossover [f|]j|. 

Figure 3 shows the kinetic frictional force, Fkin, as a 
function of the velocity v x for W = 10 and K = and 
1.1. For comparison, Fki n for W — and K = 1.1 is plot- 
ted in the inset. In the present model F^ in is given by 
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where the contribution from the viscous term in Eq. (1) is 
subtracted and the angle bracket means temporal average 
p6|| . In the low velocity regime, v x <C 10, the periodic 
potential plays no role and the velocity dependence of 
Fkin is weak. In the zero velocity limit, both FkmS agree 
with F c 's for K — > in Fig. 1. The weak velocity depen- 
dence of Fkin may relate to the fact that the impurity 
force acting on two-dimensional CDWs is independent 
from the sliding velocity. |I3] As v x increases, the veloc- 
ity dependence crosses over to the behavior, F^m oc w" 1 . 
When v x is increased further, the dependence on the 
periodic potential amplitude, K, appears. This behav- 
ior is consistent with the perturbation theory 
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and the growing order along the longitudinal direction 
observed in Fig. 2. The main contribution to the ki- 



2 



netic frictional force comes from the dynamical lattice 
deformation along the longitudinal direction, u x . (|^] As 
discussed later, the periodic potential becomes effective 
to the kinetic frictional force only when the character- 
istic length of u x grows sufficiently In the low velocity 
regime, u x fluctuates so strongly due to the random po- 
tentials. Then, the characteristic length is small and the 
effects of the periodic potential almost vanishes. When 
v x increases, QLRO of u x grows. Hence the fluctuation is 
suppressed, and then the lowest-order perturbation the- 
ory works. It is easily shown by the perturbation theory 
that the kinetic frictional force is given by the sum of the 
contribution from the periodic potential and that from 
the impurity potential. Both contributions have the same 
velocity dependence in the large velocity regime, they are 
proportional to v' 1 . Il4|jl7| -[l^| As a result, in the large 
velocity regime, the strength of the kinetic frictional force 
depends on the amplitude of the periodic potential, but 
its velocity dependence in that regime does not depend 
on it. The periodic-potential dependence of the kinetic 
frictional force is the manifestation of the crossover from 
MTG to MBG. 

We next investigate the transverse response of sliding 
lattices. Giamarchi and Le Doussal developed a kind of 
the scaling theory on the MBG state for <i-dimensional 
lattices under impurity potentials. (3|^] They clarified 
that if QLRO exists along the transverse direction, then 
a certain pinning force works in that direction even in 
sliding states. According to their analysis, the transverse 
depinning threshold force F y is given by 

in the large velocity regime, where rij is the impurity den- 
sity. For numerical evaluation of F y , transverse force F y 
is applied to a sliding lattice driven by longitudinal force 
F x , and then transverse velocity of the lattice is moni- 
tored in steady states. Figure 4 shows i™ against the 
longitudinal velocity v x for W = 10 and several values 
of K much smaller than W . Finite transverse thresh- 
old force is clearly observed. For K = 0, F y shows a 
crossover behavior. In the low v x regime, v x < 3, F y is 
almost constant. When v x is increased, F y becomes pro- 
portional to v x 2 . Such a crossover behavior is also pre- 
dicted by Giamarchi and Le Doussal around the velocity 
v x ~ F c in the present unit, which is consistent with the 
present result. In the limit of vanishing velocity, v x — > 0, 
F ' y is much smaller than the depinning threshold force in 
Fig. 1, F c . Such reduction of the transverse pinning force 
in comparison with the longitudinal one may come from 
kinematical effects. For finite K 's, i.e., in the presence of 
the periodic potential, peculiar velocity dependence ap- 
pears. For finite but small values of K , K = 0.6 and 0.7 
in Fig. 4, F y shows a similar behavior to that for K = 
in the low and intermediate velocity regimes, and seems 



to follow the scaling relation. As v x increases further, 
however, F y takes a minimum value at certain velocities 
and then increases. The velocity at which F£ is minimum 
decreases with increasing K. For larger values of K , F y 
is constant in the low velocity regime and then increases 
monotonically. In the large velocity limit, F y approaches 
the depinning threshold field F c in Fig. 1 in the absence 
of impurity, W = 0, for any values of K. This means 
that, in the transverse direction, sliding lattices in the 
large velocity limit feel only the transverse barrier by the 
periodic potential U p . 

These behaviors are explained as follows. The trans- 
verse pinning is essentially related to the transverse lat- 
tice deformation u y . df| in the case of weak periodic 
potential and the low velocity regime, the pinning effect 
mainly comes from impurities. The characteristic length 
of u v along the longitudinal direction is small. It is to 
be noted that the discommensurate structure is indis- 
pensable for the pinning by the periodic potential in in- 
commensurate systems, which are investigated here. The 
discommensurate structure of u v is responsible for the 
transverse pinning by the periodic potential. There are 
many discommensurate structures with different length 
scales for a given boundary condition. They are regarded 
as excited states from the ground-state discommensurate 
structure and the state with larger length scale has less 
excitations and larger pinning energy. When one of the 
characteristic lengths of u y along the longitudinal and 
transverse directions are much smaller than a typical 
length scale of the discommensurate structure, which will 
be the order of the neighboring discommensuration dis- 
tance, the resultant discommensurate structure has much 
excitation and a small pinning energy, and then the ef- 
fect of the periodic potential to the transverse pinning 
almost vanishes. Hence that effect vanishes in the low 
velocity regime. As the sliding velocity increases, the 
characteristic lengths of u y along both directions grow. 
The discommensurate structures with larger length scales 
can be realized there. They have a larger pinning energy, 
and then the effect on the transverse depinning thresh- 
old field is strengthened. In the large velocity limit, the 
state has the ground-state discommensurate structure in 
the absence of impurities and the transverse depinning 
threshold field approaches the depinning threshold field 
F c in that case. The velocity scale where the periodic po- 
tential becomes effective to the transverse pinning should 
depend on its amplitude, and is large for the weak po- 
tential. As a result, in the case of a small value of K, F y 
decreases with increasing v x and then increases as the pe- 
riodic potential becomes effective. When K is increased, 
that velocity scale decreases, and then the velocity at 
which F y is minimum also decreases. Further increase 
of K makes the velocity scale smaller than the crossover 
velocity of the impurity transverse pinning, and then F y 
increases monotonically with the sliding velocity. This 
scenario is consistent with all the results presented in 



3 



this paper. 

Note here that, when the impurities are absent, the 
present model reduces to two one-dimensional lattice 
models along the x- and y-axes, respectively. In the 
present study the external force is applied mainly along 
the x-axis. Some features of pinning and sliding may de- 
pend quantitatively on its sliding direction to some ex- 
tent and on the lattice symmetry. However, the essential 
features in the presence of both periodic and random 
potentials will be unchanged, even though the direction 
of external force and the lattice symmetry are changed. 
That is, in the presence of both potentials, the frustra- 
tion that causes the reduction of F c will occur generally, 
and the crossover of i™ is expected to appear for any 
anisotropic periodic potential if the underlying periodic 
potential causes certain finite static pinning force origi- 
nally in the transverse direction. 

In summary, we have investigated the static and 
dynamical behaviors of two-dimensional elastic lattices 
under periodic and impurity potentials using a two- 
dimensional Frenkel-Kontorova model with impurities. 
It is shown that the impurity can assist the depinning. 
The crossover from moving Bragg glass states to mov- 
ing transverse glass states is clearly seen in the structure 
factor. The change of the lattice order affects the kinetic 
frictional force and the transverse depinning threshold 
force. The latter shows a peculiar velocity dependence 
in the presence of both periodic and random potentials. 
Such a transverse pinning effect may be observed also in 
experiments for a vortex lattice of type-II superconduc- 
tor and vortices in a disordered Josephson junction array 
under a magnetic field. 
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FIGURE CAPTIONS 

Fig. 1. Longitudinal depinning threshold force F c ver- 
sus the strength of impurity potential W , for several 
strengths of periodic potential K 's. 

Fig. 2. Structure factors in sliding states. (&)v x = 0.96, 
(b)v x = 31.6, where W = 10 and K = 0. 

Fig. 3. Kinetic frictional force F^in versus velocity v x , 
for several strengths of periodic potential K's. 
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Fig. 4. Transverse depinning threshold force F% versus 
longitudinal velocity v x , for several strengths of periodic 
potential K's. 
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